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Abstract 

We present a relativistic constituent quark model to analyze the mass spec- 
trum and hadronic properties of radially excited u and d quark sector mesons. 
Using a simple Gaussian function as a trial wave function for the variational 
principle of a QCD motivated Hamiltonian, we obtain the mass spectrum 
consistent with the experimental data. To do the same for several observ- 
ables such as decay constants and form factors, it seems necessary to include 
both Dirac and Pauli form factors on the level of constituent quarks. Taking 
into account the quark form factors, we thus present the generalized formulas 
for the rho meson decay constant and form factors as well as the 7r7 transi- 
tion form factor. We also predict several hadronic properties for the radially 
excited states. 
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I. INTRODUCTION 



Based on a detailed analysis of the meson mass spectrum by Godfrey and Isgur (GI) [|TJ, 
the light-cone (LC) approach has been adopted to describe the pion decay constant 0, the 
charge form factor of the pion 0-0], the electromagnetic form factors of the p meson || 
and the radiative pn and urn transition form factors |3|,(|]. In this approach, Cardarelli 
et al. used the eigenfunctions of the effective qq GI Hamiltonian as orbital wave functions 
which consist of a truncated set of at least 38 harmonic oscillator (HO) basic states. While 
these wave functions reproduce the meson mass spectrum very well, it became necessary 
to introduce Dirac and Pauli form factors on the level of the constituent quarks (CQ) in 
order to reproduce hadronic meson properties. This has been done to calculate the charge 
form factor of the pion, the radiative pre and ujtt transition form factors (both including 
Dirac and Pauli CQ form factors) and the electromagnetic form factors of the p meson || 
(Dirac CQ form factor only). 

Recently J7|, however, it was observed that both the masses and the hadronic properties 
of ground state pseudoscalar and vector mesons can fairly well be reproduced by taking just 
a single lS'-state HO wave function given by 

^ 2 ) = ^7^ ex p(-2^)- ( L1 ) 

In these calculations, the choice of the Gaussian parameter (3 has been obtained by the vari- 
ational principle for the QCD motivated Hamiltonian including the Coulomb plus confining 
potential. In this paper, we adopt this model |F| and extend it to the radially excited states 
for the u and d quark sector where the masses of u and d CQ can be taken equal, m u = m^- 
The necessity of CQ Dirac and Pauli form factors was also confirmed in our analysis of both 
ground state and excited states. We thus present the generalized formulas for the rho meson 
decay constant and form factors as well as the 7T7 transition form factor taking into account 
both Dirac and Pauli CQ form factors. Following the ground state analysis, it may not be 
so unreasonable to take 25* and 35* HO wave functions given by 

= + ^ exp {~w) (L2) 



and 



2v / 30vr 3 / 4 / 3 11 /2 V ^ ^ / * \ 2p 

for the first and second excited states, respectively. However, we tried the variations from 
these wave functions and obtained more justifications for using Eqs. ( |1.2| ) and ( |1.3| ), as we 
will present in the next section. 

The LC formalism (for a recent review see ||), first introduced by Dirac |J, presents a 
natural framework to include the relativistic effects which turned out to be crucial to describe 
the low-lying mesons JTJ. Distinctive features of the LC time (r = t + z/c) quantization 



include the suppression of vacuum fluctuations [10] and the conversion of the dynamical 



problem from boost to rotation [11], which has the compact parameter space. Various 
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hadronic properties are calculated using the well established formulation in the Drell-Yan 

West q + 



frame 12 



We also use the Melosh transformation |T3[ to assign proper 
quantum numbers J PC to the mesons. 

In our model, the meson state \M) is represented by 



\m) = *2q\qq), 



;i.4) 



where Q and Q are the effectively dressed quark and antiquark. The model wave function 
consists of the radial wave function <p n s, the spin-orbit wave function 1Z obtained by the 
Melosh-transformation and a Jacobi factor 



*% Q = *(x,k ± ,X q ,X q ) 



dx 



TZ(x, k±, X q , \g)(p nS (x, k±). 



1.5) 



Here, the Jacobian of the variable transformation {k n , k±} — > {x, k±} is given by ^ 
ix{i-x) anc ^ ^he s pi n - or bit wave function for spin (S 1 , S z ) is given by 



K s ,sAx, k±,X g , Xg) = J2( X 



XX' 



lZ\f(x, k±, m q/ 



x ( — , XX'\SS Z 
,22' 1 



A A 



T^ M {1 - x, -k±,m q ) 



where the Melosh rotation is given by 



TZ t M (x,k ± ,m) 



with n = (0, 0, 1). Also, the wave function is normalized as 



m + xM — ia(h x k) 
m + xM ) 2 + k\ 



A' 



;i.6) 



;i-7) 



;i.8) 



where k = (k n , k±) and k n = (x — 1/2)M . 

We treat the meson mass according to the invariant meson mass scheme where the 

meson mass square Mq for the case m q = rriq is given by Mq = ^z^y with m q being the 
CQ mass. The intrinsic LC variables x and k± are given by x = k q / P + = 1 — k^ / P + 
and k± = k ql _ — xP± = —k q ± + (1 — x)P±. Here, the subscript J_ indicates the component 
perpendicular to the LC quantization axis n and the + component of a 4-vector k = (k°, k) 
is given by k + = k° + h ■ k. The total 4-momentum of a meson in the frame of P± = 0j_ is 
therefore given by P = (P + , M 2 /P + , 0j_) and those of the quark and antiquark are given by 



k i ) and k„ 



x)P+, 



respectively. 



' xP+ ' - 1 - / "* ' i ~ \ \ I ' xP+ ' 

The paper is organized as follows: In Sec. we fix our model parameters. First, in 
Sec. II A , we set up the QCD motivated effective Hamiltonian for the qq interaction. Using 
the wave function in Eq. ( |1 . 1| ) we determine the Gaussian parameter (3 from the variational 
principle. With a relativistic hyperfine interaction we then determine the best potential 
parameters and CQ mass m q to fit the mass spectrum for the radially excited states of the 
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u and d quark sector as well as the pion decay constant /„■. In Sec. |IIB| , we introduce the 
CQ form factors and fix the parameters associated with these form factors using the pion 
form factor F n and the transition magnetic moments fi pn and \x^. In Sec. |TI, we present 
the generalized formulas taking into account the CQ form factors and predict the rho decay 
constant f p , the rho electromagnetic form factors and the 7r° — > 7*7 transition form factor 
F^. We also predict several hadronic observables involving radially excited states in this 
section. Summary and discussion follow in Sec. TV[ In the Appendix, the details of the 
fitting procedure for the CQ structure parameters are presented. 



II. FIXING THE MODEL PARAMETERS 



A. CQ Mass and Potential Parameters 



We start with the QCD motivated effective qq Hamiltonian as given by 
description of the meson mass spectrum 



IH for the 



H, 



ss z 

nlm 



{Ho + V, 



qqJ 



ss z 

nlm 



M, 



<?<? 



nlm I 1 



(2.1) 



where M q q is the mass of the meson, the free Hamiltonian is H = ym 2 + k 2 + ym? + k 2 , 
l^nim > is the meson wave function as given in Eq. (|1.5|) and V qq is the addition of the 
central potential Vq and the hyperfine interaction Vh yp that we present in detail below. 

As shown in 0, we use the variational principle for the ground state wave function in 
Eq. ( |1 . 1|) to determine the value of the variational parameter /3 by satisfying 



d{</> ls \(H + V )\<j> ls ) 
dp 



0. 



(2.2) 



where Vq is the interacting potential consisting of (a) Coulomb plus HO, and (b) Coulomb 
plus linear potential 



V (r) =a + V con{ (r) - 



4k 
3r 



(2.3) 



Here, V con f(r) = br 2 [V CO nf( r ) = br] for HO [linear] confining potentials. Once the optimal 
parameter j3 for the ground state is determined, we use the same j3 for the first and second 
excited states. This ensures that all these states are orthogonal. 

To distinguish the vector meson from the pseudoscalar meson, we include a hyperfine 
interaction fUJll 



E q Eq 



9m q mq 




E q Eq 



(2.4) 



m 



I + k 2 , E B 



+ k 2 and (S q ■ Sg) = 1/4 [—3/4] for vector [pseudoscalar] 



where E q - 

mesons. Since we are dealing with light mesons including radially excited states, this rel- 
ativistic correction is essential. If we were to calculate the mass splitting AM n s between 
pseudoscalar and vector mesons using the non-relativistic \4 yp given by 
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«»m - S^ vVc °-' = (2 - 5) 

then one gets 

AM 1S : AM 2S : AM 3S = 8 : 12 : 15, (2.6) 



where 

AM nS = ( n5 



327TK 3 . . 

<T(r) 



9m q mq 



(2.7) 



This result clearly contradicts with the experimental values AM^ = 0.630 GeV, AM 2 $ ~ 
0.100 GeV and AM 35 « 0.100 GeV. 

To determine the potential parameters, we proceeded as follows. First, we chose the 
quark mass m q = mq as an input parameter assuming m u = ma- For both HO and linear 
potentials, we tried reasonable CQ masses in the range 0.150 GeV ^ m q < 0.300 GeV. 
The potential parameters a, b and k were chosen to provide an optimal fit for n, tt(2S), 
p and p(2S). Then, we predicted the meson masses for ir(3S) and p(3S). Our results are 
summarized in Table |. As one can see, the results for a wide range of the CQ mass m q are 
within the experimental limits of the mass spectrum. Furthermore, the difference between 
the results of HO and linear potential is quite small once the best fit parameters are chosen. 
Note that the identity of the higher resonances is not completely clear yet [IS|. While there 



are some experimental evidences for rho meson with mass 1700 ± 20 MeV, there is also 
a reported resonance at 2149 ± 17 MeV. In fact, p(1700) might be a SS'-hybrid mixture. 
Further consideration should be made to see if our simple model may or may not be suitable 
to reproduce these higher resonances. 

We have also examined the Gaussian smearing function to weaken the singularity of the 
Delta function in the hyperfine interaction 



5 3 (r)-^^expHxV). ( 2 - 8 ) 



It turns out, however, that the smearing effect is negligible in our model calculations and 
the well known value of a — 1.8 [1|] did not change our results appreciably. 

The radial wave function in Eq. (|1.1|) has been used successfully to approximate the 
ground state wave function in a couple of papers as mentioned above. In order to make sure 
if the first and second excited states can be well approximated with Eqs. ( |1.2| ) and ( |1.3|) , we 
varied the wave function of the first excited state to the mixed wave function 

|2S) = / 2 |2S)+/ 3 |3S), (2.9) 

where \2S) and \3S) denote the second and third HO states. We determined the parameters 
f 2 and fs using the variational principle in Eq. (|2.2|) as f 2 = 0.982 and f% = 0.190, showing 
that the added term is much suppressed. For the second excited state we used the wave 
function 

|3^ = / 3 |2S)-/ 2 |3S r > (2.10) 
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which is orthogonal to Eq. ( |2.9| ). However, the calculated mass eigenvalues for these states 
deviate less than 2% from our values determined by using Eqs. (1^2) and (|1.3|) . Therefore, 
we trust that our approximation of using Eqs. ( |1.2| ) and ( |1.3| ) is well justified. 

While all the parameter sets summarized in Table | give a good agreement with the mass 
spectrum, the pion decay constant is rather sensitive to the choice of m q . Following the LC 
approach of [|T6l , |T7i for the pion decay constant given by 



one gets 







97 759 



P,00)V2P+ = iP + V2f n , 



(2.11) 



(2tt) 3 / 2 



d k i dx 



ldk n 



dx 



m? q + k\ 



--(f>xs(x, k±). 



(2.12) 



0.190 GeV and the HO 



As shown in Table |, the calculated value of 91.9 MeV for 
potential is in good agreement with the experimental value of 92.4 ±0.25 MeV from [ IB] . We 
therefore from now on use the parametrization where m q = 0.190 GeV and f3 = 0.4957 GeV. 
In our work, we don't need to introduce an arbitrary axial-vector coupling constant on the 
level of CQ @. 

In Fig. |l|, we compare our central potential with those of other publications |],[14],]7] in 
the interesting range of up to 2 fm. Our potential seems quite comparable to the other 
calculations. 



B. CQ Form Factor Parameters 

While the CQ were often treated as pointlike particles, there are hints in the litera- 
ture [0-0] that the CQ may not be treated as pointlike particles. Especially, the Gerasimov 
sum-rule calculation [jTj| indicates that the CQ should be treated as extended particles. 
Thus, we introduced Dirac and Pauli form factors on the level of CQ. In fact, the Gaussian 
parameter (3 = 0.4957 GeV of our parametrization is somewhat larger than that of other 
calculations [|7|, |I1||IT| . |T4"| where (3 is generally in the range j3 ~ 0.3 . . . 0.4 GeV. Therefore, it 



seems important to take into account these CQ form factors in order to obtain comparable 
results with the experimental data for the hadronic properties. 

We substitute j M = e g 7 M at a quark-photon coupling by the more general form 

J g » = F^\Q 2 )e ql , + F { P q \Q 2 )K q ia, u £-, (2.13) 

where the Dirac and Pauli form factors on the level of CQ are normalized as F^(0) = 
Fp\o) = 1. Here, e q and K q are the CQ charge and anomalous magnetic moment, Q 2 = —q 2 
is the 4-momentum transfer square and a^ v = |[7^, j u ]- For the Dirac and Pauli CQ form 
factors Fjj and Fp, we adopt simple monopole and dipole forms 

F^\Q 2 ) = -L and F P «\Q 2 ) = - \ _ (2.14) 

l+<rg } >QV6 (\+<rf >Q*/12) 
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as other authors j| used. It turns out to be sufficient for our purpose to use this simple 
monopole and dipole form. 

To fix the model parameters < r 2 D >, < rp >, k u and Kj, we fit our model to the available 
pion form factor data and to the experimental values of the transition magnetic moments 
fi p7T and /i W7r . The details of the fitting procedure are presented in the Appendix. 

Our determined values for the parameters < r D >, < rp >, k u and are shown in 
Table O- Note that these parameters are comparable to those determined in H. Also, we 
calculated the ratio (e u + « u )/( e d + K d) which is predicted in ]18| as —1.80 ± 0.02, which is 
comparable to our value of —1.94. 

In Fig. 0, we show our calculation for F n and compare it to the experimental data taken 
from [ 19[| . We also compare our calculation to a simple Vector Meson Dominance (VMD) 



model where F^ MT) (Q 2 ) = 1/(1 + Q 2 /M 2 ). In Fig. 0, we show the radiative transition form 
factors F p7T and F W7r and the body form factors Hp and H P n . 



III. CALCULATION OF HADRONIC PROPERTIES 

Having fixed all the parameters of our model in Sec. |J, we now present our predictions 
of various hadronic properties of pseudoscalar and vector mesons including also the radially 
excited states. 



A. p° Decay Constant f p 



From the definition [17[ 



<0\J qii \P,lJ 3 > V2F< 



(3.1) 



we calculate the rho decay constant f p including both Dirac and Pauli form factors on the 



level of CQ by taking J q ^ as defined in Eq. fl2.13j ). We obtain for p° 



where 



[uu 



dd)/V2 



f P =^ {[e u - e d }F { D q) (M 2 )I D + [k u - K d ]F P q \M 2 )I I 



(3.2) 



and 



V3 



(2tt) 3 / 2 
A(27i) 3 / 2 m q 



d k\ dx 



x{l 



d k\ dx 



m + 



2k 2 



4k 2 



A 



x(l — x) \ A 



>is 



'x,k\ 



(3.3) 



M <f>is{x,k ± ). 



(3.4) 



M 2 
p 



Note here that the momentum transfer square equals the mass of the rho meson Q z 
(0.770 GeV) 2 . While Id has been already derived in [17], Ip is our new body form factor 
related to the Pauli form factor in the electromagnetic current operator J q ^. Using our 
parametrization (Table |I|) we get the value f p = 153 MeV which is in a good agreement 
with the experimental data f p = 152.8±3.6 MeV obtained from the width T(p — > e + e~ 
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B. p + Form Factors 



Our calculation for the p + form factors follows the one presented in [20|. However, we 
again treat the CQ as extended objects characterized by Dirac and Pauli form factors in 
contrast to a treatment as point like ones. 

In the standard LC frame the charge, magnetic and quadrupole form factors of a meson 



can be obtained from the plus component of three helicity matrix elements [21 

1 



F c 
Fm 
Fq 



(2a + 1) 
2 

(2a + 1) 
1 

(2a + 1) 



1 4^-\(2a-3)F^ + l(2a-l)F; 
o y 2a o o 



(2a - V 



2a 



F + - 



F + - F' 



a + 1 



a 



(3.5) 
(3.6) 
(3.7) 



where Fy X = (P', X' |(J£ — J%)\ P, A), a = ^2 is a kinematic factor and is defined 

in Eq. ( |2.13| ). This representation is not unique; there are different prescriptions in the 
literature, as discussed in 0. This ambiguity is reflected in the fact that the angular 
condition 



A(Q 2 ) = (1 + 2a)F^ 



F^_ - V8aF+ - F+ = 



(3.8) 



is in general violated unless the exact Poincare covariant current operator beyond one-body 
sector is used. 

At zero momentum transfer, these form factors are proportional to the meson charge e, 
magnetic moment /ii and quadrupole moment Q\\ 



F c (0) 



eF M (0) = 2M pf i 1 and eF Q (0) = M 2 Q X . 



(3.9) 



In the LC quark model, the matrix element {P , A Jg P, Xj can be calculated by the 
convolution of initial and final LC wave function of a meson 



P\ X' 



PA 



E 



d k i dx 



' dk n / dk 



dx V dx 

U(kq>,Xq') J qU(k q ,Xq) 



-Mx, k ± )<f>ls(x, k' x 



lh+ 2 

K q , 



kn 



Tl'*(X q/ , A ? ; 1, X')K(X q , A 5 ; 1, A). (3.10) 



After a straightforward calculation, choosing the + component of the current, we get 

F^(Q 2 ) = (c - e d )F^(Q 2 )I x D \Q 2 ) + (k u - K d )F ( P q) (Q 2 )lP(Q 2 ), (3.11) 
where the body form factors related to the Dirac part are given by 

1 



i°J(Q 2 ) 



d k\ i dx 



Idkn dK 
dx V dx 



(j>is(x,k±)(j)l s (x,k' ± 



x(l — x)AA' 



x 



1 - 2x) 2 k ± k' ± + (2x[l - x]M' Q + m) (2x[l - x]M + m) 



(3.12) 



S 



d 2 k 



d.r\l ^ L \l^ IL (j)is(x, k±)(j)l s (x, k' ± 



1 - 2x 



x 



2xk x M' (M' -M ) + Q (2k 2 y - 2x(l - x)M' M Q - M' m 



(3.13) 



i + D -{Q 2 



d 2 k 



<l r\l^\l^^(j) is (x,k ± )(j)* is (x, k' ± 



X 



dx V &r ™ ' ' xy M M^(l-x)AA / 

(ml„ + mMg + 2m 2 + 2x(l - x)M M^) k x k' x + (2m 2 - 2x[l - x]M M^ k. 



-2k x k x 



m 



(X'kl + Xk' x 2 ) + 2k* 



(3.14) 



and 



I + d + {Q 2 ) = I d 2 k ± J dxs^^MxMW ls {xA) 



M M^( 1 -x)XX' 
(k'\ + mX') (k 2 ± + mX) + (k' ± k ± f - (1 - x) 2 Q 2 k 2 y 

+ (2m 2 + mM + mM' Q + [2x 2 - 2x + l]M Mg) fc^jfej 



(3.15) 



Here, A' = 2m + M' Q . These formulas have already been used by the authors of || where 
they calculated the p + form factors including only the Dirac part of the CQ form factors 
by assuming that the anomalous magnetic moments of the CQ are negligible. We, however, 
present a complete calculation including also the Pauli form factor part. For the new body 
form factors related to the Pauli part, we obtain 



/°°(Q 2 ) = / d 2 k ± I dxJ^J^Mx,k ± W ls (x,k' ± )- t Q[] 2r) 



dx V dx 



2x(l — x)mXX' 



x 



mX' + 2k' ± 
Mi 



~k x 



mX + 2k\ 



M 



k .... 



(3.16) 



d 2 k ± I dx y~^V~^^s(x, k±)(j)* is (x, k' ± ) 



Q 



2\/2x(l - x)mXX' 



x 



1 — 2r / -> - 

— ([m + xM^]k'±k± - [m+ (1 - x)M' ][k' x k x - k 2 } 



M' n 



mW ( [Am + 2kl][K2 ~ k " ] + [x ' m + k ' ±2][Xm + 2k 



(3.17) 



d 2 kt 



dx\ 



I dk„ / dk 



dx V dx 



Q 



X 



(m + [l-x]M^) [k' x 
(m + [1 - x]M ) (k 



/2 



2x(l - x)mXX'M M' Q 
2k x k 2 v ) + (k' ± 2 + A'm) (m + xM )k x 

k\ + Am) (m + xM' Q ) k' x ] (3. 18) 



„2 

2k x k 2 
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and 



d 2 k 



dx\ 



'dk n dki 



<Pis(x,ki-)<P*is(. x ^' 1 _ 



Q 



X 



dx V dx ^^•"^^•"^2x(l 
-(m + (1 - x)Mq) (k\ + Am) k' x -(m + xM ) (k x 



+ (m + (1 - x)M ) (k'j 2 + Am) k x + (m + xM' ) (k' x 



i m 
k x ~~ k y ) + 2k x k y 



+ 2k x k 



(3.19) 



We show the result for the p + form factors using our parametrization (Table ||) in Fig. |j. 



In Table |n, we list our calculations for /i X and Q\. Note that our result is quite comparable 



to other calculations [22,6,20 



Using Eq. ( |3.8| ), we have also calculated the violation of the angular condition and 
compared it in Fig. |5] with several other calculations |6l,|20"|j . The violation of the angular 
condition seems to be suppressed by taking into account the structure of CQ, compared with 
the other calculations. 



C. 7T° — > 7*7 Transition Form Factor 



The 7T — > 7*7 form factor F nl has been calculated quite successfully in several models 
where the CQ has been treated as a pointlike particle jl7|,|7|. We show the generalized 
formulation including both Dirac and Pauli form factor of the CQ. 

The 7T° — > 7*7 transition form factor in leading order is defined as 



iG q {Q 2 )^ pu P v e?q*°, 



(3.20) 



where P is the momentum of the incident pion and q* is the momentum of the virtual 
photon. Taking P + = 1 the vertex factor is given by 



V3 



J d 2 k± J dx\dx25(l — x\ — X2)J -7^-4>is(x,k 



v I A,A',A 

x 7e(A,A;0,0) 
1 



dx 

vx(x 2 , -k±) ux'(xt, k± + qj_) ux>(xi, k± + q±) + ux(xi, k ± ) 



j: 



Jl _ (fc ± +g x ) 2 +"t 2 k\+m? 
9-L H 



{xi <-> x 2 ) 



■>-■> 



(3.21) 



Note that we have to plug in J q ^ as defined in ( |2.13| ) at the vertices of the real and the 
virtual photon. Here, of course, the momentum transfer square for the real photon should 
be zero. After a straightforward calculation, we get 



F^{Q 2 ) = Gu{Q 2 ) - G d (Q 2 



(3.22) 



where 
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Gg(Q 2 



V2Q{2tt) 3 / 2 



d k i dx 



X 



M 



' M' 2 x 2 {l -x 

K, 



-4>is(x, kj 



e q F D (Q z )— mQe g - (k x + 2(1 - x)Q) xM 2 + (1 - x)k x Q< 



+n q F P (Q 2 



mil 

2j 



+ejl - x]Q fc^e 



x)M 

Q r 



2m 

e q Q xk^ + m 2 +e g [l— x]Q 
kl + m 2 -(l-x) 2 Q 2 



k y 



(3.23) 



Our result is shown in Fig. || As Jaus [17] pointed out, the hadron structure of the neutral 
pseudoscalar meson (jr°) may not be well enough approximated by the one-loop calculation. 
Gluon-exchange effects may introduce additional structure that might lead to a mechanism 
analogous to the flavor mixing of isoscalar states. Especially, at Q 2 = 0, from the one-loop 
formula given by Eqs. ( |3.22|) and ( |3.23| ) and the definition of T ni given by 

7T 



rr ~ - ja 2 F^fM, : 



3 



(3.24) 



we obtaion V 



7T~f 



4.97 eV. However, the agreement with the experimental data IV 



7.8 ± 0.5 eV can be obtained by taking into account the PCAC and the anomaly of the 
axial- vector current |23]24], which predicts 

^(0) = 7Z^- ( 3 - 25 ) 



Using this and Eq. ( |3.24j ), we obtain r ?r7 = 7.82 eV. 

On the other hand, the high-momentum transfer region (Q 2 > few GeV 2 ) is dominated 
by the one-loop off-shell quark contribution as evidenced from the Q 2 -behaviour (~ ^2^. 
Recently p5| , it has been shown that for high momentum transfer region the transition 
form factor can be obtained by the renormalization scale and scheme independent PCQD 
calculation as 

2U (, 5 ay(e- 3 / 2 Q) \ 

Q 2 { 3 7T ) ' 



^3.26) 



where ay(e 3 ^ 2 Q)/n fa 0.12. Using this relation, one can get a good agreement with the 
experimental data for Q 2 > 2 GeV 2 , as shown in Fig. 0. 



D. Predictions for Radially Excited States 

We also calculate several hadronic properties of radially excited states using the same 
parameters shown in Table ||. For these calculations, we use the radial wave functions <p2S 
and 03s from Eqs. (|1.2| ) and ( |1.3| ), respectively. 

In Table [IV], we summarized our results including the transition magnetic moments n p >+ n + 
and /Vti-o, which are the transitions from the first exited state to the ground state. While 
the pion and rho decay constants are proportional to the wave functions at the origin for 
the ground state, the nodal structures of 25 and 35 radially excited states make intuitive 
predictions on the ordering of decay constant magnitudes untenable. We also show our 
result for the pion form factors for the first and second excited states in Fig. |7|. 
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IV. SUMMARY AND DISCUSSION 



In this paper, we used a simple relativistic CQ model to calculate various properties of 
radially excited u and d quark sector mesons. We find that a QCD motivated Hamiltonian 
with a relativistically corrected hyperfine interaction yield the mass spectrum comparable 
with data. We also treated the CQ not as pointlike objects but as extended ones, thereby 
introducing both Dirac and Pauli electromagnetic form factors for the CQ, and were able 
to obtain a reasonable agreement with the data for hadronic properties like the pion decay 
constant /„., the pion form factor F^ and the transition magnetic moments /i p7r and /i W7r . 
Furthermore, we used this model to calculate the rho decay constant f p , the rho form 
factors and the 7r° — ► 7*7 transition form factor F nl , presenting for the first time generalized 
formulas including both Dirac and Pauli form factors on the level of the CQ. We find that the 
angular condition A(Q 2 ) is better satisfied when the CQ form factors are included. However, 
our model is intrinsically limited to small momentum transfer range. As an application of 
our model, we also calculated several hadronic properties for radially excited states. Further 
experimental data on radially excited states will give more stringent test of our model. 
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APPENDIX: DETAILS OF THE DETERMINATION OF THE CQ PARAMETERS 



The calculation of F n has been considered in a number of references [26,27[fl[4||. For 

7r + = ud one gets 

FAQ 2 ) = (euF^(Q 2 ) + e^(Q 2 )) H^Q 2 ) + (k u F p u) (Q 2 ) + k s F^(Q 2 )) H p (Q 2 ) 

= F%\Q 2 )Hl + (k u - K d ) F p «\q 2 )H p (Q 2 ), (Al) 



assuming the same internal structure for the u and d CQ, Fp = F^' = F^' and Fp 
F p d ^ = F P q \ The body form factors H D and Hp are given by 0] 



HKQ 2 ) = J d 2 k ± J dx\ 



ldk n dk' n t r m 2 + k ± k' ± 



and 



(A2) 



H P (Q 2 " 



d 2 k 1 / dx 



dkn -l d % ls (x,k ± )<f>l s (x,k' ± h ^ 



dx V dx 



2xMM 



(A3) 



Mq is given by M' 2 = ^gf and k' ± by k' ± = k ± + (1 - x)q ± . 
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Following a similar spin and flavour algebra for p + = ud, u = l/y2{uu + dd) and 
7T° = 1/v / 2(mm — dd), one gets for the radiative p + — > 7r + 7* and uj — > 7T°7* transition form 
factors 



F P .(g 2 ) = (eui#>(Q 2 ) - erf(g 2 )) H P J{Q 2 ) + (k u F p u \q 2 ) - K d F®{Q 2 )) H?(Q 2 



F&HgiQ 2 ) + (k u + K d ) F¥>(Q 2 )H p /(Q 2 ) 



(<?), 



(A4) 



and 



iUQ 2 ) = (e u F^\Q 2 ) + e.-F^g 2 )) ^(Q 2 ) + (k u F p u \q 2 ) + k 3 F?\Q 2 
= F$H%{Q 2 ) + {k u - K d ) F^(Q 2 )H^(Q 2 ), 
where the body form factors Hp and Hp n are given by [ffl 



2 / d 2 k i / rfx 



dx V 9a; 



is 



, 7 mX + 2k 2 y 
x, k , 21 



(A5) 



(A6) 



and 

^r(g 2 ) 



d & i / da; 



'<9A;„ /<9fc 



<9x V dx 



-cj) lS (x,k ± )(j)l s (x,k , ± 



X (m 2 + k±k' ± ) - 2M (1 - x)k 



x(l - x)M M{ ) mX 



(A7) 



Here, X = 2m + M . 

We are now able to fix the parameters k u and K d . The value of the transition form 
factors at Q 2 = 0, the transition magnetic moments /i pvr and )U W7r , have been experimentally 
determined from the radiative decay widths of p and mesons p!7|j20|] viz. 

r2\3 



T(p -> vr 7 ) = iaF p ^(0) 2 



2M„ 



and a similar expression for V(u — > 717) as p p7r = 0.741 ± 0.038 GeV 1 and //, 
0.06 GeV -1 , respectively [15|]. From (|A4|) and ( |A5|) we get 



and 



(A8) 
— 2.33 i 



^ = -^i + («„ + Kd )HF(0) 



(A9) 



(A10) 



By fitting the experimental value for p p7r and p W7r one gets k u — K d = 0.138 ± 0.015 and 
Ku + Kd = 0.036 ± 0.01 or, correspondingly, k u = 0.087 ± 0.013 and n d = -0.051 ± 0.013. 



The mean square radii (MSR) associated with F^ and F P q> are < >= — 6 



dp. 



(?) 



dQ 2 



and < r*- 9 '' 



2 

p >= 



-6 



dF 



(?) 



Q 2 =0 



dQ 2 



. Assuming the same internal electromagnetic structure for 



(«) 



quark and antiquark, we have r D 



r" D ' = r^ and r p 



(a) 



(?) 



rp' = r P 
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To determine and r y p' , we fit our model to the well known electromagnetic form 
factor for the pion F n , following the procedure outlined in H. Using ( |A1| ), we get 



< 



dQ 2 



< r D > + < r H * > +(k u - K d ) < r H n > 



(All) 



where < rlU >= 



_5 dH p 



dQ 2 



and < r 2 >= -6 g 



are the MSR associated with 



the body form factors Hp and ifj. Using our values for k u and Kd, we find < r 

0.432 ± 0.016 fm 2 



0.0495 ± 0.022 fm 2 to obtain the experimental value < r 2 >- 
last parameter to be determined, 



28 . The 



is not affected by the Q 2 ~ range of F n (see 
Eq. ( |A11| )) and can be chosen to be < r 2 p >= 0.136 fm 2 to get a good agreement in the 



< rf, >, 



whole range of existing pion form factor data. 

Finally, making use of the uncertainties due to the experimental errors, we determine 
our parametrization (see Table |l|) to get the best possible fit for the rho decay constant f p . 
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FIGURES 

FIG. 1. The meson interaction potential Vo(r). Our parameters are compared with the HO 



potential of Jj]], the quasi-relativistic potential of ISGW2 (k = 0.3 and k = 0.6) [14| and the 
relativistic potential of GI Q. 

FIG. 2. The 7r + form factor times Q 2 . We also show the predictions from a simple VMD model. 
The experimental data is taken from [19| 

FIG. 3. The radiative p + — ► ir + j* and uj — > ir°j* transition form factors F pn and F^n- We also 
show the body form factors Hp and Hj? 

FIG. 4. The rho form factors including both Dirac and Pauli form factors on the level of CQ. 

FIG. 5. The angular condition tested by A(Q 2 ). Our result is compared to other calculations 
[20,|]. A feature of our model is the fact that the violation of the angular condition is smaller than 



in other models. 

FIG. 6. The it — > 7*7 form factor i ? 7r7 . Our result agrees with the experimental data pi^pOfl 
only for low Q 2 . Using a simple power law Ansatz as suggested in [25] we get agreement for 
Q 2 > 2 GeV 2 . 

FIG. 7. The ir + form factor times Q 2 for the ground state (IS) and the first two radially excited 
states (2S and 3S). 
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TABLES 



TABLE I. Model parameters for the best fit of the mass spectrum for both HO and linear 
potential for several quark masses m q . The parameter j3 has been obtained from the variational 
principle for the ground state. For comparison, the parametrization from together with the 
ground state masses are shown. Also shown is the pion decay constant f n used to distinguish 
between the different parametrizations. All units are in GeV, except the parameter b which is in 
GeV 3 [GeV 2 ] for HO [linear] potential, f n which is in MeV and k which is dimensionless. 





HO 


HO 


HO 


HO 


HO 


HO 


linear 


linear 


HO a 


linear 3. 


Experiment 13 


m q 


0.150 


0.180 


0.190 


0.200 


0.250 


0.300 


0.150 


0.250 


0.250 


0.220 




a 


0.268 


0.264 


0.256 


0.246 


0.243 


0.220 


0.112 


0.187 


-0.144 


-0.724 




b 


0.0127 


0.0117 


0.0118 


0.012 


0.0092 


0.0069 


0.089 


0.05 


0.010 


0.18 






1.207 


1.22 


1.22 


1.218 


1.246 


1.26 


1.165 


1.242 


0.607 


0.313 




P 


0.4782 


0.4894 


0.4957 


0.5019 


0.5224 


0.5463 


0.4914 


0.5277 


0.3194 


0.3659 






0.136 


0.139 


0.138 


0.138 


0.138 


0.137 


0.139 


0.138 


0.137 


0.135 


0.135. ..0.139 


M <2S ) 


1.280 


1.272 


1.275 


1.279 


1.261 


1.253 


1.288 


1.258 






1.3±0.1 




1.799 


1.797 


1.806 


1.817 


1.802 


1.805 


1.781 


1.776 






1.801±0.013 


M p 


0.770 


0.770 


0.0770 


0.770 


0.770 


0.770 


0.770 


0.770 


0.770 


0.770 


0.770±0.0008 


M p(2S ) 


1.437 


1.441 


1.448 


1.456 


1.459 


1.470 


1.442 


1.457 






1.465±0.025 


M p(3S ) 


2.051 


2.049 


2.059 


2.071 


2.061 


2.069 


2.033 


2.037 






c 


fn 


75.8 


87.9 


91.9 


95.8 


113.9 


129.9 


76.2 


113.9 


92.4 


91.8 


92.4±0.25 



Predictions from [7]. 
b Data taken from [15|. 



"See the text for a discussion of possible resonances at 1.700 ± 0.020 GeV and 2.149 ± 0.017 GeV. 



TABLE II. The model parameters. 



m q = 0.190 GeV < r 2 D >= 0.050 fm 2 

a = 0.256 GeV < rf> >= 0.136 fm 2 

b = 0.0118 GeV 3 k u = 0.074 

K = 1.22 K d = -0.048 
P = 0.4957 GeV 
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TABLE III. Observables for different model parameters. The first column shows the calcula- 
tions from [^t]]?]]. The second and third columns show our predictions without and with CQ form 
factors. 



Choi/Ji Hg,@| our w/o CQ FF our w/ CQ FF Experiment 



m [GeV] 


0.250 


0.190 


0.190 




/? [GeV] 


0.3194 


0.4957 


0.4957 




<r 2 D > [fm 2 ] 








0.050 




<r 2 p > [fm 2 ] 








0.136 




Ku 








0.074 












-0.048 




U[MeV] 


92.4 


91.9 


91.9 


92.4±0.025 


< r 2 > [fm 2 ] 


0.448 


0.327 


0.427 


0.432±0.016 


^Anomiir -> 7*7) [eV] 


7.73 


7.82 


7.82 


7.25±0.23 


f P [MeV] 


151.9 


231.7 


153.0 


152.9±3.6 


M^+lGeV- 1 ] 


0.782 


0.603 


0.701 


0.741±0.038 


^^[GeV- 1 ] 


2.35 


1.81 


2.27 


2.33±0.06 


P : Mi 


2.2 


2.27 


2.63 


a 


p ■ Ql 


0.20 


0.53 


0.69 


a 



a There are no data available yet. However, the results of other theoretical calculations are given 
by: /xi : 2.26, 2.1, 2.3 and Qi : 0.37, 0.41, 0.45 from B, |(J and fH, respectively. 



TABLE IV. Shows predictions of several hadronic properties for radially excited states. p p i+^+ 
and jU^/^o are transition magnetic moments for the transition from the first excited state (25) to 
the ground state (IS). 
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25 


35 


U [MeV] 


91.9 




18.6 


42.3 


< r 2 > [fm 2 ] 


0.427 




0.606 


0.854 


f P [MeV] 


153.0 




90.5 


96.4 


A^+IGeV- 1 ] 
/WGeV- 1 ] 


0.70 




0.55 


0.48 


2.27 




1.79 


1.58 


/y+^IGeV- 1 ] 




-0.23 






/VTroIGeV" 1 ] 




-0.71 
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